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MoO LRACT 


This theory permits the calculation of pressure 
CiemmeGemEmons Over 2 tiin airioil in steady, inviscid, 
incompressible flow; or given a desired chordwise pressure 
distribution the camber line of an ideal wing can be 
determined. Ihe method treats the circulation about the 
wing :-as a continuous vortex sheet of variable strength 
covering the wing planform and trailing downstream to infinity. 
itemec Mee sOsmeocePrecsene numerical solution for the pressure 
distribution solution are not satisfactory. However the 
solution for the camber line of an ideal wing are in: reasonable 


agreement with published two dimensional results. 
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eR ODUCTLON 


This analysis is a continuation of two dimensional 
thin airfoil theory [ Ref. il and lifting line theory [Ref. 2| 
into a three dimensional lifting surface theory, using the 
concept of a continuous vortex sheet over the wing which 
trails off to infinity in the downstream direction. The 
continuous vortex sheet differentiates this method from 
the vortex lattice kernel methods of Multhopp [Ref . 3 | ; 
Falkner [Ref . 4 | , Weissinger [Ref. 5 | , and Watkins et. al. 

[Ref. 6 | and double-lattice finite element methods of 
Lopez and Shen | Ref. 7| and Giesing et. al. [Ref . 3 | ; 
The present method is restricted to steady, inviscid, 
incompressible flow and thin airfoils which are symmetrical 
in semispan. 

BeeCmemectlac On abeul the Wing is treated as generated 
By a continuous vortex sheet of variable strength covering 
the wing planform and trailing downstream to infinity. The 
vortex sheet strength is set by the Kutta condition at the 
trailing edge and requiring no flow through the wing at 
specified control points. 

The flow equations for the wing are resolved into down- 
wash angles, with assumed series solutions for the circulation 
functions. From the known planform geometry and the required 
solution at the control points, no flow through the wing, the 
constants of the circulation functions series are determined. 
The equations are put into matrix form, the matrices being 
formed by numerical summation over the planform, and solved 
Miceh@arceebra. Once the constants of the circulation 
Peer noiomare Known all of the airfoil characteristics can 
be calculated. 

This analysis generates the spanwise and chordwise 
pressure distributions over the wing which are necessary in 


designing the wing structure. The solution is also done in 


10 





Pe emcepeENadte to the pressure distribution over the wing is 
specified and the shape of the mean camber line of an ideal 
airfoil is determined. 

This analysis uses the consepts of additional and 
reference lift which 1s of academic interest in under- 
standing airfoil theory. 


Lol 





Il. General Problem Development 


Peete VORTEX SYSTEM 

The present analysis uses the theory of thin wings to 
describe the characteristics of lifting surfaces in steady, 
inviscid, incompressible flow. The most common planforms 
are those which involve straight leading and trailing edges, 
as illustrated in Fig. 1, but the present analysis applies 
also to wings having curved leading and trailing edges. We 
consider only wings that are symmetrical with respect to the 
x axis. We also restrict the present analysis to lift 
distributions which are symmetrical with respect to the x 
AX1S o 

The wing and the trailing vortex system associated with it 
are treated as a continuous vortex sheet of variable strength 
ane sheet strength ¥ is regarded as a vector in the 
X,y plane. Ihe wing vortex sheet is treated as lying in or 
very close to the x,y plane and a 1s everywhere tangent 
to the vortex sheet. 

The sheet strength ¥ is related to the velocities Vu 
and Vw just above and just below the vortex sheet at a 
given point x,y and to the remote velocity Ve , as indicated 


in the vector diagram, Fig. 2. From the vector diagram Fig. 
Ae 


Vo = Vet $e (1.1) 
Vi = Ve-sY¥%k (1.2) 


It then follows that 


—= = —» 


il yu + = ae Gere) 


(We JV.) = Yxk (1.4) 


WZ 





or conversely 


=e ———in. 


we kx ( ye Va GPRS 


In Fig. 3, consider the circulation dJ’ around contour 
AOBB O AA. Line segment AOB lies just above the vortex 


sheet, and segment BO’ A lies just below the vortex sheet. 





— 


A0oB= AO 8B = ds | (1.6) 
Then = 

ire § Vide = \ ERs [ae [te 

P= Veds +O -VUeds +0-(U-V) eds 

df=(¥xk)-ds (1.7) 


It is assumed, which will subsequently be confirmed, that 
Maem an equation (1 7 ) is an exact differential. There- 
mere, there exists a function 


‘ie xn) 


such that for any small changes in the coordinates 
or of 
d= (Se-)dx + (Sy) dy a 
Translating this into vector terms gives: 


gp (73 +722). (Td + Tdy)= vk (1.9) 


Comparing equations (1.7) and (1.9) shows 
yi! Y x le Chale 
Or conversely that 


Y= kx VC Giraiel) 


ES 
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Xr 


Wing Coordinates 
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baler 
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Fig. 2 Velocity Vectors 





Fig. 3 Vortex Sheet Element 





These results confirm that dl! is an exact differential. 
Translating equation (1.11) into Cartesian components 





gives 
b= - (25) au 
(er (se) Ce sy) 


The vector function A defined by equation (1.11) or by 


scalar equations (1.12) and (1.13) is non-divergent. 


o” 
=a) 





y & omy _ Oo (— 
See 


- ole 
a 5 wr )+(2-)= O Coles, 


OY 


4 
1S) ty od al 


Newton's Law applied to a fluid element of fixed mass 
PAK a y42 1S 


= = eax ayar BY OF) 


{it is assumed that the fluid is non viscous; that is, shearing 
stresses are absent. The resultant pressure force of the 
fluid element arising from the static pressure variation in 
‘the fluid, is the negative of the pressure gradient multiplied 


by the increment volume. 
resultant = ( - grad p ) Ax44 42 


The pressure force plus the weight is the external force, 


therefore the equation of equilibrium becomes 


i DV 
_ grad = ae (272) 
Og = G P= aa | 
which is Euler's equation. 

If the fluid is assumed incompressible and equation (2.2) 
1S integrated dropping the unsteady term it becomes 
Bernoulli's equation 


via 
Cer +P- C42 = constant (203) 
16 





The pressure difference across the wing at an arbitrary 


point can be found from Bernoulli's equation. 


(PnP). (AP). ()* (EY 


aver VV.) . 2 Ie Rk) 
Via - vo 
_. 25. (i Sees ee at) 
: od ‘ Vow , Veo \ . 
AP). ¥y 
ce = 4 (+4) C2) 


Equation (2.4) ultimately fixes the lift distribution over the 
wing surface once the circulation function J’ has been found. 
In the trailing vortex region behind the wing, there is 

NoOlwine Surrace, and hence no pressure difference across the 
vortex sheet. According to the Kutta condition, the pressure 
difference drops to zero not only behind the wing but also all 
along the trailing edge. This is expressed as follows: 





& ¥ 
(42). 2 2 (2 Jy = 4+(s#)_=0 al 


This result shows that the vortex lines behind the wing and 
along the trailing edge are all parallel to the xX axis, with 
no spanwise components. Spanwise vorticity [ Y4 can occur 
only along the wing surface. 

In Fig. 4, next page, consider the incremental lift dt 
exerted by the doubly shaded portion of the wing, that lies 
between y and (y+tdy) and between x and x. Thus: 


Xx 
tf “ 
ris be Ve | (SE) dx dy 
Ru 


i 





i 


: 
dh Veo | [ 25 dx | dy 
dL = 0 Vel dy (2.6) 


This result is a form of the Kutta-Joukowski Law which states 


ni 














a 0.9%, OO I IIL IIL LLL LAS Ld 






Pewee Lite Lategration 


ED 


that the force experienced Dy a body in a uniform stream 1s 
equal to the product of the fluid density, stream velocity, 
and circulation and has a direction perpendicular to the 
stream velocity. Here [™ represents the circulation about a 
GOmeGmeme im tae Plane y = constant. Contour C passes in 
HEOQMemenmEine Wile and Plerces the wing Suriace at point x,y. 

If we now extend the upper limit of integration to the 
trailing x-, moldtmag 7 COnStant, we Obtain: 


na 
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X—> Xr 

(x4) — TF (¥) 

eee 

Ji, = CAS Bag ORD) 


thus é Gp etixes tae spanwise lift distribution over the wing. 
In order to reduce equation (2./) to a convenient dimension- 


less form, it is useful to introduce the following auxiliary 
notation. 


2 





S = total wing area, ft C28), 
Pose = wine span, ft (2.9) 
a = = mean chord, ft "(2 505) 
me bk. : (2.11) 
a. aspect ratio sii 
Ge = at = spanwise coordinate OBE), 
dg= 33 (2.13) 
dt ne 
Cy = -EcoCCEronelitmuicoeniieient (2.14) 
Gat du 
C= <= = wing lift coefficient CAS) 


Pate liziie tae above nomemclature, it 1s easy to reduce the 
spanwise lift distribution as given in equation (2.7) to 
either of the following two equivalent forms. 





(Es d= ba (S)dy= (se) en (2.16) 


To find the total lift we integrate over the span, For 
Sees real Iift distribution it 18 only necessary to 


integrate over the semi-span and multiply by two. Thus: 
1S, 





f 


Cre (a () de = 2{ (= )da (2.17) 


° © 





See DUCED VELOCITY 

The velocity dwp induced at an arbitrary point P, 
coordinates Xp on the wing surface by an element of the 
ieeeeeoleet wiich dies at the point whose coordinates are 
x,y is found as follows. tJIhe vortex sheet element has length 
ds, width dn and strength ee inemerreularion strength of 
the vortex filament involved is Jr = Ydn. 


The Biot Savart Law for this case may be written in the 


form: 
de> Teta ds (3.1) 
where 
meee eg 4) + & (2-2, ) (3.2) 
dn ds = dS = element of the wing area 
VU}. ee Geen Val (3.3) 


Neglecting the terms &(Z-Ze) in equation (3.2) and RYz in 
eamation (3,3). 





= eS ar fh Ih 

y= - Coe, ssa (Ca (368) 

also _ a . 

L \ fe 

rxy = |(x-xe) (4-4) OO 

ok Nie 
5 52 0 
~ Se 

= ke [(x-%)(2E) + (4> Ye) (Sq) | (3.5) 
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Fig. 5 Induced velocity in the x-y plane 


21k 





Upon substituting expressions (3.2), (3.3), and (3.5) into 
(3.1), and integrating over the entire wing surface S$ 

( but not yet over the vortex sheet trailing behind the wing) 
the downwash due to the wing vortex is obtained. 


x-Xe) (s=— el 4) (Sp s | 
Wp = pe | ee = =) a 


HS 
(x-% - 2 (yesae ar (3.6) 


The incremental downwash velocity dw, induced at point 
P due to the semi-infinite trailing vortex filament of strength 
dl> which leaves the wing at the point whose coordinates are 
XT and y as shown in Titers by the Biot Savart Law: 





d[+=- (22) Jy 
dw ~ ee [ | |- Ci Xe) 
aa | ae ) [ (x-xe)** (4 -4e)? re “ 


Integrating this over the entire span of the wing gives: 





i, ig (X-%p ) Ae 
= ————_s ——— sy tear egeenicnanete Al ; 
_ =m | Fol Ffonnse Cogs] | Tidy G8) 


The total induced velocity at point P is the sum 


We = We + We (3.9) 
where Wwe and We are evaluated respectively from equations 
Peo sndeee es) eine first of these intesrals represents the 
velocity induced by the vorticity over- the actual wing surface, 
and the second integral represents the velocity induced by 
the vorticity trailing behind the wing. 

The downwash integrals are now non-dimensionalized and 
the Cartesian coordinates x and y are transformed into 
dimensionless relative coordinates 3 and @ defined as 
follows: 


Lee. 








CN 
iH 
iH 

NH 
cy 
to 
at) 
Gr 
ft’: 
<a 
OD 
Q 
ey 
O 
H 
fat 
= 
f- 
” 
@ 
@) 
‘e) 
© 
by 
O 
ieee 
| 
Qs 
al 
4) 


C 10) 


4 | 
(@e— 9 = Belative spamwise coordinate Shisle, 


C= = = average chord (Sul) 
A = = = aspect ratlLo r18) 
SS fy) eS , , 

AY = = leading edge function (rie 

x) = = = trailing edge function Cr > ) 
a = ( (7) = relative chord function (372.6) 


Git Seana) - a7) (3.17) 


Note that if the parameter A, is specified and if any two 
Sumeie three fumctions involved in equation (3.17) are 
specified, this information completely defines the wing 
planform. 

The circulation function f’ is non-dimensionalized as 
follows. Ihe number four is inserted into these definitions, 
Somewhat arbitrarily, because this reduces the final downwash 
integrals to a particularly convenient numerical form. 


a Pr , , , 

P ( 3) = ries = dimensionless circulation 
function along the wing 
surface (3 eke.) 

® 

l (7) = mio = dimensionless circulation 


4 Vos function at and behind the 


trailing edge (3.19) 


ae 





° ° . am / —- e « e 
The induced velocities We and Ww, are non-dimensionalized 
to corresponding downwash angles x,- and e&,s as shown below. 
Those angles are defined as positive if down from free stream 








velocity. 
’ We -k . . 
“Ap = V = downwash angle point P induced 
i Ey VOrELeClLeyeGiStri bution Over the 
wing CS 5 2OL) 
uw Wee lo . * 
Re = — 7 = downwash angle at point P induced 
ia DyvGmEleity distri pution) behind 
the wing Crs 
Ke = Fee, = total downwash angle induced at point 


P by the complete vortex sheet system 
of the wing (B22 ) 


Some useful geometrical relationships are 


(x-xe) = 8 HL (K_+C2) — (Xie + Ce 3p) | (3.23) 
(ul=Glaee 2 (eee (3.24) 


Also at any point on the wing 


meee ves F (35 7) (3.25) 
Demterentiating equation (3.25) 


® 


P 
eg. + (25) ay - 4 Vos [ ($5) d i+ (37) del (3.26) 





Differentiating equations (3.24) and (3.25) 
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GaZz7) 


Mee s =! Cdz+(K. +C z)gn-0] 
(3.28) 


dy= sd 


The prime marks in equation (3.27) denotes differentiation 
Substituting equations (3.27) and (3.28) 


with respect to y. 
into equation (3.26) gives 


- 2 [Cd +(Ki+Caldr]: 


sg 
ee a, = (3.29) 


oe dy 


SE 
: ig (Sa) da 


The coefficients of dz and dy must be separately equal on 
Poumectdeawor equation (3.29) since this relation must be 
. Hence 


Soelomeeamnor arbritrary values of dz and dy 





















(-) ZC = te | af (3.30) 
P , Min ai 
(Sr) (EL + C2) + (Sy) + te Gq) 9) 
Solving those two equations for a and 2 gives 
Bo) = te ae (SF) 7) 
Hui C4) (30) (3.33) 


Se ee Siw 
Sa) = tM [ (se) ( € 
The element of wing area in relative coordinates may be 


expressed in the following manner 
25) = Ix), duJB = 09] aa 
ME) degq=S Bclszeq G3) 





dS= 2s° 


Z2 





ce 


Fig. 6 Element of Wing Area 


The Various quantities of equations (3.6) and (3.8) are 
mEamsmommead Linco non-dimensional notation. 
5/2, 


[ (ame) elamsey] = 9°] ge [(5er CD-R Cb] 0-7 I" 


[(x-xe) 2S + (4-40 )so]gS= Ves S aL [(X,+Cf) - (Keele) 
-)> ere) (2 2) — ( Sith) (20) ts ZCd5d4 





& (S3 


= $V +f 2] (Ze Ci) -(Ke+ 6 5-)] | | 
[ 2~ Ze) )[e BF) (XU C4 (SE) ] f 2 


+ 2| 
FR 
= Ht ZEl( K+ Cf)-(Zee Coe] -(-e) Ke CH] GF | 
- (2-20) C(3E) g ¢4 Cree 
Bee Ce 7- ) | (3.37) 


(Xr -%)=s & & |x ~ | 
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Study of the algebraic structure of these relations 
shows that the final downwash integrals can be reduced to a 


relatively concise form by adopting the following auxiliary 
nomenclature. Let 


2{ (K.+C4) ~ e+ Ce ge) -(1-Xe) (Xu+ C% \] 





eee cee Ne (3.38) 
(a [Be Cf) - Se *Cefe) | + (4-Re) : 
F, = z(-7e) C ioe 


| , 
Ge [xr- (Zin eCrie) |e PE | CP — .40) 


The downwash angle integrals equations (3.6) and (3.8) can 
now be expressed in the following form. 


1 I 
“ 


Xe ar | {FE titer (3F1 +Fp (LU be EN eh de (3541) 


-| Con 





ke 
_ Let d ft 
Ke = at F. 5 tole) Oe )dn (cua) 
D. THE CONDITION OF NO FLOW THROUGH THE WING 


The shape of the mean wing surface near the x,y plane 
may be expressed by the function 


z= z(x,y) CD) 


A unit vector mre normal to the above surface at the arbitrary 
point P on that surface, in terms of its Cartesian components, 
is given by the expression 





x ~(se)ec 7 Sk +(\k 


- | (2 » S1GSeqeh sal ale Cae) 


; pz \2 OZ 2 ee ; . ; 
Since (Se and aoe are negligible in comparison with’ 





unity, equation (4.2) is adequately approximated by the 
linearized version 





= an re oZz\ = aa 
a erat oe SSS) Se i Cap) 


The resultant velocity vector at point P may be closely 
approximated by the expression 


Vo= © Vg COSK + KV SINK — BR We (4.4) 


Equation (4.4) implies that the induced velocity we is taken 
as positive if down. 

Since x is a small. angle equation (4.4) may be linearized 
CO 


eee (a =e) (4.5) 
The condition of no flow through the wing now requires that 

as =O (4.6) 
thus 

fc Vet + E ( oo -We)| Eve: eo) i i S=\+ RW -9 (a7 ) 
OAS 

— Vo (S=\, + (Voo% -we) = O (4.8) 


hence 
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oz = _ W/ 

ep > ava (4.9) 
sjbee 

WwW, 7 ‘ ry, 

. an (4.10) 
therefore 

2 f ’ 
(2S), = X- (Ke + x; ) (4.11) 


The overall wing equation is now found by substituting 
the integrals of equations (3.41) and (3.42) into equation 
(4.11). The result is the basic wing equation. 





) - Bae | | fr. (SF) FRE) ds dy 





¥ + fests 25) )dy iD) 


tievumrivence functions F,, Pos and F, are as defined by 
equations (3.38), (3.39) and (3.40). These three functions 
depend only on the planform and aspect ratio of the wing. 
lijteemucrrtediiy, equation (4.12) must be satisfied at all points 
in the wing surface for an exact solution. However, an 
adequate approximate solution can frequently be obtained 

by satisfying equation Ce) ac a sutiieient number Of 
discrete control points suitably distributed over the wing 
surface. 


Pee NCE AND ADDITIONAL LIFI DISTRIBUTIONS 


Let the wing slope function (22). be known at every 
point 3, , 4, Of the wing. Also let the entire leading 


ES, 





edge of the wing lie exactly in the x-y plane and the 
Paileenes edee lle in the x-y plane at midspan. In general, 
however, Z2 is not necessarily zero at other points along 

the trailing edge. Let eas 7) be a function which 
Selerpes equation (4.12) when x has a particular value ar. 
The subscript r stands for reference condition. Thus at the 
reference condition equation (4.13) becomes 


, 

SE 

Sele =e tea f [ [PGE 
-{ 0 








)+Fe(2#) dian 





t x 

+e { Fs (SE) dy (5.1) 
=| 

The general solution ore or equation (4.13) can now be 

written as the sum of the reference solution and an additional 

Sserueron in the form 


P'Uge) = 03,4) + wK-me) fo. (3.7) (5.2) 


To demonstrate that equation (5.2) is valid, and to obtain 
the equation that governs the form of the additional solution 
function we proceed as follows. Substitute equation (5.2) 
feo equation (4.12). Substraet equation (5.1) from the 
result. Finally divide through by the factor (a«-a,), The 
result is 


Fs 


sete te ff (hE) GR) fers 








te c, (Site) Jy Gap 


Thus equation (5.1) governs the reference solution or the 
particular solution while equation (5.3) governs the 
Pagmuetonal solution. Ihe specific value of Ay in equation 


ee) ts fixed by the particular condition that 
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oe? = ee 7 = 0 (5.4) 
Pamven Chiat 
Z2—O0at 7 = 0 tor ail 7 Sis), 


Ihe reference solution depends not only on the wing 
Peomrorm and aspect ratio but also on the wing slope function 
and the angle of attack aw, . However, the additional lift 
function involves only the wing planform and aspect ratio. 

It will be shown later that the wing lift curve slope 
@epends only on the additional lift function [a and is 
independent of the reference solution [. . Therefore, 
the wing lift curve slope depends on wing planform and 
aspect ratio only, and is independent of the wing slope 
Pele t1 ON . 

To specify the design of a wing it is necessary to fix 
two planform functions x4) and cC(7) , the aspect ratio 
AW , and the wing slope function ——- ene higi erealjes 
pe, Ye on the wing. For a wing of specified design, equation 
(5.1) must be solved for the initially unknown function 

im 3 7) and equation (5.3) must be solved for the initially 
Geknown function [x (3,77) : 

In the case of equation (5.1) the above procedure can 
also be partially reversed, that is, the function [7 (j,7\ 
can be specified arbitrarily, within certain limits, and the 
equation cam then be solved for (23). as the unknown. fhis 
reversed solution procedure had two advantages. Firstly, it 
means that the lift distribution over the wing in the 
memeEcMeescondlt1Om Can be stipulated independently and that 
the wing slope function can always be found such as to yield 
the desired reference lift distribution. Secondly, it 
Pelg@ems that equation (5.1) is much easier to solve when the 
wing slope function is the unknown than when the reference 
lift distribution is the unknown because, in the latter case, 
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the unknown function is under the integral sign. If (22), 

is the unknown, equation (5.1) can be solved by direct 
integration; usually this must be performed numerically. 
Repeated integrations of equation (5.1) are required, a 
Separate numerical integration being needed for each discrete 
point at which the wing slope function is being evaluated. 

In order to represent the wing slope function adequately 

it 1s necessary to choose a sufficient number of discrete 
points suitably distributed over the wing. 

Unfortunately, the above reversed solution procedure 
cannot be applied to the additional lift distribution ae 
as governed by equation (5.3). Of course, the equation can 
still be solved, at least approximately. 

The choice of the reference angle of attack &r which 
fixes the corresponding reference solution Na Gls ce 6) eeezulia 
extent arbitrary, The most common choice is to set ex, equal 
to the angle of zero lift of the wing, that is to the angle 
miamen ylelds a wing lift coefficient of zero. For a 
symmetrical lift distribution the lift coefficient is 


CLs ( & (=) dn oa 
o : 

The restriction that the 1ift coefficient equal zero does 
not necessarily mean that the integrand is identically zero 
all along the span, just that the net area inder the curve 
sums to zero. Thus at Cr = Q, portions of the span may be 
exerting upward forces, while other portions are exerting 
downward forces. A wing that behaves in this way at Cr = Q 
may be said to be aerodynamically twisted. Conversely a 
wing that at C,; = 0 yields Cy-= identically zero for all 
values of 7 may be said to be aerodynamically untwisted, 
this however does not imply the absence of camber. 

The reference condition for this analysis 1s that which 
corresponds to the angle of attack for which the wing lift 
Coerficient equals zero. tJlhis particular lift distribution 
is also commonly labeled as the basic lift distribution. 
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F. ANGULAR COORDINATES 
a further shift from the linear relative coordinates 
2, to corresponding angular coordinates @ ,¢e , is now 


made. See Fig. 7. The basic conversion relations are 








3 =% (1 - coss) (6.1) 
4% = cose Cony) 
then 
S = % sin ¢ ens) 
— = Sant CS (6.4) 
ok )= xB, A (6.5) 
Oe cae — 





Fig. / Angular Coordinates 
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The downwash integrals can now be converted as follows 





T & 
— Xp = ae {f n arr ey )- Fe ciel So) £m ¢ CSHedg do 
mam © 

















7 re 0 
af dlr o I+ 
nae af Glee $f QB de oo 
© ° 
where 
G, = Fy sin 8 (6.9) 
Gy = ae sin d (6.10) 
Ga = = F. Goan) 


Now by refering to the above relationship and to the previous 
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equations which define the functions By Fos and P3 Le is 


possible to summarize the final wing equations as the 


culmination of a series of systematic calculation steps. 


The following scheme of auxiliary notation is used. 


Depending on planform, the derivatives 


Cee 


Ne 
ie 


hip 


h 


we 


= cos ©& 


= cos 9P 


= x<,(e)+ C(s) 3 


= A, ( Sp) + C (ee) Ap 


* ( 1 - cos ¢@) 


(1 - cos %& ) 


- & [x.o]- 2 ge [cel] 


dX de 
de ad 4d © 


may be multi-valued at midspan, however for a symmetrical 


planform they are set to Zero. 
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hs 
Cas = Sl aehic 














0 ~ %e (6.26) 
The basic wing equations become the following 
Reference lift 
nN ww 
SZ 
Sele = Xe + tam | [Lo RE) ] déde 
+ J os (SE) Altes (6.27) 
Additional lift 
{f H 
Die a1o 
o=t+te | ( [6 (Se)+ G2 SSll ded 
Le) en) 
I aes 
~ £ [Gs (S t=) Je CGieziey 


Gs mPFROXIMATE SERIES SOLUTION 
Ihe reference and additional distributions can be adequately 


represented by double series of the following form. 


RH OS 
r. (go, e) Paes H-(¢) SINGS Cg 
To £21,3,S-- 
Ro oS 
fa(d,e)= > aes H-(6) SiN Se (7.2) 
C=Q SA,3,S 


For a spanwise symmetrical lift distribution, only odd values 
of s are used. Ref. 1. There are (R41) distinct functions Hr 
and ( ==: ) distinct sine functions in the above series. 


The total number of initially undetermined constants in each 
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of the above series is therefore 


=(R+1)( = ) GES) 





If there are (R+1) chordwise stations and ( =<. ) spanwise 
stations over the semi-span, this establishes T control 
points. By satisfying each governing equation at each of the 
T control points we can solve for each set of T constants. 
Once the constants Brs and Ars are found, all aspects of the 
wing performance can be readily calculated including 11ft 
distributions, slope of the lift curve, induced drag, etc. 
The character of the H functions is established as 
follows 





a= & (22) (7.4) 
oa 2 FPR Vow ary (7.5) 
d* C: o% 
therefore 
- (46 e- Ae oa (+8 | -= SING (RE C76») 


For the additional lift, equation (7.2) is substituted into 
Bamation (7.6) 


RH Ss 
ao GC <5) ) ) Ars sae (gir )[swse Gham 
Ho S=t3,5- 


For the reference lift, equation (7.1) is substituted into 
equation (7.6) 


Rel oS 
A Sie [AS sf) SINS © 
a a cd Sans me ag) | (7.8) 


It is now possible to make use of the results of thin 
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airfoil theory and set up the fimctions in equations (7.7) 
and (7.8) in an analogous form. Ref. 2. Thus for the special 
Sacewor ~- = 0 we use 


SHo_ 2 dHo = 1+ Ug 
SF =6SNG dF = samy toes 


and for r # 0 we use 





oHr = Ze alr 
Equations (7.9) and (7.10) can be integrated to yield 
Eee =O 
Ho ($)= 4 (p+ SND) (7.11) 


Bom c= | 


Hi (d)= (gp — S22) (7.12) 


fom r= 2 


He (6) = fj] Sees _ see Cre) 


Veral + | 


At the trailing edge, G =tr, these relations yield 


Uy in = = (7.14) 
fe) * (7 15) 
Pein = © (7.16) 


By substituting these values into equations (7.1) and (7.2) 
the corresponding expressions for the trailing vorticity 


are obtained. 


7 S 
F-(e)= E> (2Bo<* Bis\ SIN S©& (7.17) 


Sees 
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n S 
Ea (oe) > = > (ZAos + Ais) SIN S96 ise 
S5).5.5 65. 


An ideal wing is one for which (22), = (£2). =| (8) Zlie 

the leading edge. This requires that Bos = 0 for all s. 

If the planform function is arbitrarily prescribed the wing 
should not be presumed to be an ideal wing and all terms in 
Bos should be included in the solution. The additional 
lift is always characterized by the singularity (SS) + = 
at the leading edge so that the terms in Aos must always 


be retained. 
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ITI DESCRIPTION OF NUMERICAL PROBLEM SOLUTION 





peeeEC ULATION SCHEME 

Ihe calculation scheme in this analysis uses two distinct 
meshes. There is a coarse mesh of control points, and a 
fine mesh for the purposes of accurate integration. It is 
desirable with respect to both meshes that the ratio of the 
number of spanwise stations over the semispan to the number 
of chordwise stations should be roughly equal to the aspect 
ratio divided by two. This gives roughly equal resolution 
in the spanwise and chordwise directions. 

The total angular interval, 0 to Win each case is 
subdivided into equal sub-intervals and the calculation point 
located at the center of each sub-interval. The calculation 
points associated with the fine mesh must not coincide with 
any of the control points associated with the coarse mesh 
because this would cause the functions Fy and Po to assume 
an indeterminate form. While this indeterminacy can be 
resolved it is best avoided. A good way to avoid this 
difficulty is to make the fine mesh an even submultiple of 
that of the coarse mesh. Then the control points will 
always lie along intersections of the boundaries of the area 
elements which comprise the fine mesh. See Fig. 8. 

The double integration over the wing area reduces in this 
scheme to a summation over the points of the fine mesh. The 
number of terms in the series being equal to the number of 
Gemerol points. 

If the pressure distribution is to be found for a given 
Pune tne Goerflicients of the reference and additional lift 
series are solved for first and then the pressure distribution 
can be calculated. 

If an ideal wing camber is to be defined from a given 
chordwise pressure distribution, the wing equation is solved 
ior (2=), » The terms in We are dropped which eliminates 
the leading edge singularity, and an equal number of terms is 
added on to the end of the series. In other words, Bos is 
Seprequal to zero. Ihe induced drag of the wing can be 
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reduced to an absolute minimum if the terms in sin se 

are eliminated retaining. only the term in sine . This 
defines an elliptical spanwise load distribution. These are 
desirable features to incorporate into the design of an 


ideal wing, and it also simplifies the numerical solution. 


B. MATRIX FORMAT 


By analogy with equations (7.7) and (7.8) the reference 
and additional lift distributions may be written in matrix 
momma C 


ge les - [o| (3\ (8.1) 


((2).f= [ol] [A (8.2) 


The matrix format of the basic wing equations (6.27) and (6.28) 
become 


3 -xe f re ‘J | (8.3) 
fife -(R1 faq (3.4) 


The matricies ime] and [R] in the above equations are known 
T by T matricies whose elements depend only on wing planform 
and aspect ratio. 

Solving equations (8.3) and (8.4) for the initially unknown 


constants B and A gives by matrix inversion 
fe{=(r] {2 -xe| (8.5) 
-| 
Ale [a YEU ae 
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Equation (8.6) shows that the A constants depend only on 
wing planform and aspect ratio while the B constants in 
equation (8.5) also depend on the wing slope function. 

The points used in the Q matrix need not be the same as 
Exe control points used in the R matrix, but may be any 
point on the wing. In this analysis the points for the Q 
matrix were picked as six stations across the semispan with 
eighteen chordwise stations at each spanwise station. This 
permits a good representation of the functions [" and (2) & 

If the reference pressure distribution over the wing is 
specified, coefficients of the lift distribution series can 
be found, by equation (8.7). Then the wing equation (8.8) 
can be solved by direct numerical integration. 


((42).] = Lol {a (8.7) 


la > Be surg SINE SIN © (8.8) 
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es 
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+ G2> ore He cose| + az 2. Gs Bese de 
Le) 
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eo US 


Pee OULVING FOR THE WING PRESSURE DISTRIBUTION 

The numerical results of the pressure distribution solution 
Memesmoe Satisfactory. Ihe pressure distributions and circu- 
lation functions calculated by the computer program were not 
smooth or in accordance with [Ref. 1] and | Ref. 10] 5 be oe 
Pempevea this is a result of difficulties with the first 
chordwise term, Bos. the term by term downwash is plotted in 
Fig. 9 , at a spanwise station near midspan for a wing of 
elliptical planform and an aspect ratio of 20 using 5 chord- 
wise stations and 10 stations over the span. The integration 


mesh was four times finer than the control point mesh. Ihe 
first term should yeild uniform downwash. the rest of the terms 


are of the general nature expected. 
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Fig. 9 Chordwise Downwash 
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B. SOLVING FOR THE WING CAMBER 

Ihe numerical results of calculating the camber of an ideal 
wing given the desired chordwise pressure distribution were 
satisfactory. Fig. 10 compares the results of 22 -x, fOr .c: 
Wing of aspect ratio 20 and rectangular planform using a 15 
mermm rit distribution series at 5 chordwise points near mid- 
span, with the two dimensional results of [Ref. 10.| ees 
felt the differences at the leading edge can be resolved by 
fees More terms in the lift distributions Series and more 
integration points. The depicted pressure distribution was 
imposed at all spanwise stations. The results are normalized 
memo lift coefficient of one. Table | shows the results at 


Weerous sDanwise stations. 





Table 1 
aAR=20 15 term lift distribution series rectangular plantorm 
0.99 Ome9 Ore! OS OS 

9 eee RB Boe 
On025 Obs OSA: Orne oo ON sie Orly oZ 0.1817 
on206 0.0292 OS GS SIL 0.0789 OF Olere Ss! Of OILS 
e500 -0,0448 -0.0781 -0.0965 -0.1038 -O.1611 
0.794 -0.0967 ere ey S95 LS yl ear alsys -0.2297 
0.976 -0.0998 -0.1456 -0.1884 -0.2142 -0.2256 
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V. CONCLUSIONS 


The numerical implementation of the theory needs 
additional work to achieve the correct computer solutions 
for the pressure distribution solution. The areas of 
greatest interest are the assumed series solutions for 
the additional and reference lift, especially the first 
chordwise term, Bos. The camber line solution is consistent 
with two dimensional theory near midspan of a high aspect 
ratio wing and allows the camber line to be determined at 
various spanwise stations for any aspect ratio. It is 
felt that this theory has merit and deserves further study 
to correct the numerical solution of the pressure 


eieetir) DUCLON. 
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APPENDIX A. DERIVATION OF SUBROUTINES 


PemoubROULINES FOR THE PRESSURE DISTRIBUTION SOLUTION 


QMAT calcualtes the Q matrix of equations (8.1) and (8 

















“2 
AP 4 Rs y 
(<=). = CO! - ) Ars =| s25 ste | aI = (7.7) 
ae s > mia z a | SINSS (A.1) 
=}. = yas s #8/ (2, d¢ 
7 Za d Hr A.? 
HHL = Per ee “ ) 
S( = SINSS Cary) 
aie) = TA HH (A.4) 
RMAT calcualtes the R matrix of equations (8.3) and (8.4) 
(BE), cake = ote [[ (or BE +62 BE | dade 
., i 
+ {os (Sz) de (6.27) 
® RH S Ss | 
25 = >>) Bre sg ( He) Sin) SS Cy 
x YS 
SIE (A.6) 
PY, Bre Ye 
‘fic ~ d (2 Bos + Bis) SCos SS (A.7) 
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ae [or Se(Hd smise 
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oC 7 H- S$ cos st = 4) (ZBes +B] )2 Gry snise do (A. 8 ) 


E 2 - (Ane 9) 
pe ae Hr) 
Si=2 SINSO (A.10) 
Suze Vee (A.11) 
eS COs SS (A.12) 





Ger) — > > see <= |G HH1i-St + Ga- Hi. <2] 
-s CRS 25 Ae (A.13) 


Pee oUPROULINES FOR THE CAMBER LINE SOLUTION 
QMAT calculates the Q matrix of equation (8.1) which has 
Deen simplified with a single term in theta, an elliptical 


spanwise lift distribution, and the first chordwise term is 








dropped. 
eee | = 4h ys Ley cet SSS CAE) 
Fx nn eC Sw d ¢ 
ae £2. d He CSe oe 
HH = sme d¢ 
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er) = 2 - HHi- St (\.16) 


CAMBER is a numerical integration of equation (6.27) 
simplified for an ideal wing. 





ce Bg R 
SZ _ ddéde [ SINTPSIND cig 
ox OTR i a Gi 2 Br Z 


= Le (A.17) 
tT G> 2 Se H- cose a aq S GB Cas © de ° 
9 


49 





C 
C 
C 
C 
C 
C 
C 
C 
g 
G 
C 
9 
c 
€ 
C 
g 
G 
G 
e 
c 


CCCCCCUCCCCCCC COC CC ECC CUCU CECE CCC CECE CCOCCECCCCCCCC CCC CEC eC CC Ceres 


OOO OVUYU OV UOO ODO VOO0U O 


WW 
a «6D Pi 
——- = < 
— — as -_— 
Ym WWM — iL + 
mm o(/) —) Oo MO 
(a a od VY) - 
LU > “s —™ WwW AN 
el © Se —_ * 
OD <—_ ee op ee 
Iwas YW 2 & 
‘od Tg) ~i! <f 
Ie W Ze = 
Se. WW Lu 
OuWae >Z oO 
@ <f wo 2 m 
LU <f ww mm = LL fal, 
et ict Wwe Yt & 
em > Maw x 
a) mw wh nN 
CNS Rel Ww © 
ee Oe GQ) CZ x 
OER he Waa 1H Dw ‘oll 
zazCcre =e Ww MN 
Wma Tt ~2 FT iu 
m>eFeOYFr YN = —&—- &— 
ah Ww FM 
OewOo 2ew Wo 
ese) tat: => ~ 
E-=- SWE ZF NK 
Sse UE SS ae 
+ jy LL LL a QL wee wow 
YN WN waVvZay net 
—~ OL pong LL > 
_— MAWL>O li-~x 
Zeer OY Oe ONE OO 
ItweE “= ve /) VY Cc) es] 
YWOSS eOBmuiy mt je (7) 
OuUmWO waza ow a~ & 
O22 oe ZO Te ON Sm 
CIE QOQrtae—— CO 
os le | los ep Pe ee 
LQO0Or “i ™ Yam OO —— 
YW SOPethLOoOrte Fr em 
ce le a od GO) GD OO a So 
TWw—wa eo OO We eww, 
Ee> GCWertworot few Yewm 
—OeUIWWoOWT WYN srvOI 
Wow WamMtwO - 
ay aS Se * o, sa ad Pall Re FA Pe 
iat OS) = E-OOoOoO0O 
Ss —_— & WZ NMANMNY 
Seo 2 stot) Tu Fad PS GF GPT, 
Aaz OMrtDe= Wt OWWWwWW 
CL me DL ee pee LL =< ST Ts 
pa En) PP med tt et ee 
es Om i il ee oe 
TOW W rr Mew 


WwW 2VenNaanWor = 
TiL——s QawINeE- 
EFOZzOateCnETrOstO 


CCC COECCCCCC CCC EC CCC CCE CECE CECCCCCCCCCECCCECCCCE CEC ECCCC CC CCE CCE CEC CECCE 


ITB2(12)/12*0/ 


— ~ 

ay a 
INOOnNN™N 
COatk*k # OOS 
emQAIA) «© 8 « 
CAAADOO 
1 MNS 3 4 
I aem—DO OD 


—UIN AANA AT ee INO OODDO~— & 
ema NNN EE ATER QUEER A 


On ee hee en rete ee toe Lee 


El,RIB1(5)) 
E2,RTB2(5)) 


@( QSOS) HAR (AG 255 
(25) 


am, oe 
am LY) SS 
LAVAS << 
~~ eC) 
COLANY 
aN 
—_—— 
(a ee 
KWo iw 
mM @AN 


LL 


oe) 
aim, Oe oun 
CD ~~ 
ot mt 
ee es ee 
OO 0= 
~OdtS 
Tqaa<d 
amped © Lf tae LO 
ij - * ea 
& ame O™ om, 
—OOD 
CO ad mt et 
et ee ee? See 
we LOY LEV OOF 
moat ft 
KAAS 
vo) @ 2 
—_— oe ecl 
OAH 
~OD ea 
CL net end ome 
Ine =~ CE 
pm PP nee 
LU) <I me 
6. OM) 
Kroc 
-- oS 
a an SS 
Oa © <f 
ae et) 
wat ee et 
QL OD aoe 
— I) <lOm 
TAGAAoO 
a YC O~ 
~ 2 ed ww 
wOnnDEae 
Ona ETO 
mw NNO <I 
OUd eO<t 
QO AQ am =. 
xCODRS 
a © eet CO <l 
mem el 


2FOO tte LT 


% AL AED DD OD we we ELA LOA ee et SOS 


xo 


~ 
=| = pal 
onl olan) 
- — @ 
200 
-~ & @ 
OO 
on w@m 
NNO 
- @- & 
OOo 
-P- @ 
ond eed oot 
~~ & 
= al oo 
~~ = 
eat okee) 


0.0 sO 


~ = wo 


OOO 


Low YK ROA Tee we eww 
Figo oe Aen, sn I BL ee ee 


NOQOOreOF OF UO 


rr~n 


AS Pee tS St ZIFF ZZ CZZCZAA MINN 
me HE COR ee OOOO OUOOCOeamam 
ml et et re et Pet beet tet et fe pee fee be 
OW WW PPP OKT OKT OO IMN NN NNNDN NN et et te 


my Or 


= OOOK KH UH DH OK OSSSAZSZSASASHS « 


ed aed me a So Ss SD Se > Gane > i dd 
Sree CCK we YOO OQAOOO0O04l 09000 


30 








DATA XIT/1.051.0,1.0,1.0,1.0,1.0/ 


INPUT DATA 


OOO 


= 

al 

as i 
Y = 
LJ > 
= = 
~ —_ 
= —_ 
(ao << 
CO mA > 
— Wu x< 
—_ J JJ aa 
ae OL Pam pam be MN 
ae ~ 
QO — AN 
Saeed o ~.00 
~ ©§O0O0O0 aaace 
UV ew em ee 
we we LV UVLO SE we 

eae? wee geet CV) od 

qaaceaa ala} 
ada aNaead 
BSLSSLSS SRE etm mee 
Cert v LNN—we 


PRINT INPUT DATA 


WOU 


1+RP1,MESH 


p 
R 
(DZDXMACI),1=1,T) 


CCMPUTE THE Q MATRIX 


W<t 


1y,RP1,AR,Q,PHIP,THETAP) 
(CQCT,J)eJ 


CCMPUTE THE R MATRIX 


Si 


1,v),1=1,T) 


— 
= 


1,RP1,»MESH,AR oR) 
=1,T),1T=1,T) 


(( ROT, I) oJ 


SOLVE FOR THE A VECTOR 


RMely Tots AyIOGT »WK,ITER) 


(ACT) y1f=1,T) 





IDGT 


WRITE(6,500) 


a 
2) 
_ 
CO) 
UJ 
> 
oa) 
Lh — 
ole _ 
ome < 
> 
Co ~x< 
Fan 
ue OU NY 
met 
LL 
> Oem 
| — 
Oo O~ 
wy Oca 
as) 
> | 1 OD) 


SOLVE FOR THE ADDITIONAL PRESSURE DISTRIBUTION 


(CPACI),I=1,T) 
SOLVE FOR THE REFERENCE PRESSURE DISTRIBUTION 


Q:Ay,CPA,T,T;3)) 


( 
) 
) 


OOO 


Q,B,RPCyToT11) 
(RPC CI),1=1,T) 


RD ( 
O7) 
Q1) 


D2 


SP1D02+t1))*AR/4.O*PI** 2 


D2) )*AR/4.0*PI *%2 


uy 
an] 
Ida <a eS 
Sage _ 
I<d<dta< 
ce ed eo ed Cn) error Gy. Os 
CWOOUuId © 23 
ae ee QOOowud& 
~~ ORNS 
Fen 1 FON AN — 
;syb!tetm O testis { 
ALAA =O BALL uN 
HHeHHH © tt HH 4 x 
ee ee ee aa 
Q. wee wee wa tee te fee? Ee eee Nee? uw 
waite CT a Se > Xs GS b— _) 
eOOAAg F AAA Ao -O. 
aAOUOOOCO Jeat avrg =o 
HW HUW Wl ey nn WN tf wOoWH Hu 


rm Domed Pome fm fret OY et et et he 
MANO PINE OM QUINLAN EEE es 
aaidddtwew OVOOUOZ 2 YO 
OeaaactUuOato00NOO0a 
QOOQCOCVCOMACE Yr CeOOUOe 





ACAI — 
| 
LALLA 
exw HH 


Se ee | 


t atl 


WOOO 

0COa,0 

(a Ao AT a Ao 

i | oe 
a ae om ) 
en Lameed came Lemon? 
wee Saget “eee eet booed 
AIO) SPUN 
QOOVE 
O20 © 
eo, Ce Ce CD 


“oO 
om) 


SOLVE FOR THE GAMMA DISTRIBUTION 


WOUOO 


(CDS ING(K-2)*PHIP (I) )/CFLOAT(K)-2.0))-(D0SINCK*PHIP(S))/FLO 


id 

~— 

Oo. 

<— 

= 

th) — 

a6 © 

—_ e 

% Al 

a ~ 

— eae 

Ze say 

ie a 

YW — 

Se 7 

t+ -_ 

~= — a= 

Gl) ae 

“SN = 

— e ~ 

- FTN ~ 

=i] fal — ==) 

—_— —_—- ~- 

+ Ze -t 

Ob =— 4 = 

Gren). = 

=—tO | < 

ST — (5 
N Ym + 
a iaA~> ~~ > 
ate Oot NICE we OO << 
QAacda «Area = as) 
Nay botnet New i 
- oan & eA FT Qo. C= 
Sat ett) OL we i ae 
OU bE He meer LUN SE LN Lyte LULL 


Se YN PUN TN RD DD 
~~ ~~ (AN NN Nw Sree ZZ 
OO Somat yY +O em eee Od ee 
PTS il —OFrOCFrostre il Zee 
= a, it Ul bwSeeaZzZ 
QoqdtoorekwvOeocrocsatow 
| can J fam cl Ja wa Sn Joa wy is J emo > | ia | Can | |, 
od 
AND Mm SF ON 
rr Ss T+ 


HO 
tt 


55 


PlOZ+(L4+1)/2) )*DSINCL*¥THETAP (1 ))*P1/4.0 


Fe oe 

woe md 

<{ ~~ 

+h 

=~ JI 

AIS 

yy A eA, pny ~~ SS 
famed fog Fd fomd peg om <f 
a & ee & woe 
mt NO SPUN + + 
ae eet See See Se? —S}cae 
Ief<teqtet Ww we <L 
pete Ge) ~ = 
SST SS Ez BONES 
OQdqaodeq © e ata 
(od Gay Gn) <5 Gada) eg (an 9 Tew 1 Se, 
— i ib ob a mi} - eff 
ed te ey ey gay onl em, met | oom 


1] mnt pat treet tet peg Lf tees ff ee ms LL} 
Pad et et ee Set ee) a er |] ee 
FN PUINT ee ikea 
Ode de che aeaace 
TLS SSSeTSrSTSse 
eee ae Oe 
Oaqada qocnaodatu 
SIs soa OOoOO 


OQ ~ 
+ +r 


O-DCGSCUPHIPIT) )) 
(XTC 1),1= 


1 


1,18) 


PEG) SRE SUITS 





= a a — — ce7) 

mem mm MH MH OM faa) 

ee) a a ~~ = — co faa) faa) foe) fae) ~_ 
co om a cas m wn I) Oi OG ayac & ace ee me 
= - — ee & ace om a coe ~— es es ew OM 
(ad Ye «ce ec bh br) oer) WD) aD fea) 
~ ~ - ~ - - 2 O42 led aa 4 ~om a @m © amd j— 
ot =) Com — =i — & & tod & o— — -_ _— _ — 
25) ) Cae a a ao a N NN A j= —¢ tet ~ 
— Ee -_- —- - —_— dm arc) ec Oe ~ & ~ -_ 
— === Loe Hm —_ — bh iy rT AO Mm <{ 
~ oe Co od & — = (—¢ — mt ef <I <I <f <f >. 
iy Sw fas) cv} AJ amd & ~ & & 2 = = = = = 
<I x_<datiaq< <Io- (oo 2 Ie A= == fo) = <{ 
(al 3 ey 2 NR 2 1 [jee a oO = oO © Oat <I I <at<_ es) 
O @ Siaee). 6. laa: eee ice ao Ga a) a 
~ ~ - ~ ~ ~ Zw oe w tw we a = mg « — 
— mtg — » » ~ ~ 4 r—| mettle = 
< x K KK x = = — =o x K KX x ~< 
~ & o~ a ~ x“ >< x< < < « od a a ~~ ~ 
= =! ol onan = nd ~ ~ ~ ba ment aan] — om = —= 
oe a ey fh es oO -— —- ~— - - 6. oO 7 Oo. a; fo 
a rea « x a ~aeenete « 
LS ~as — ete —w C) — we ew — — re = 
o al ‘el jen {ee GF se Ce el fel fal (oWek To a G. G. 
OEZNHZNSZMOMSHPSEMOMMOSVNSNAHSNEMOSSTVNSBVNSHHSuTEMSsErsZsEMmn=z 
ae<tbe gk << B < <eop <ok S ece e ee  e eeOp eg Oe 


NO SNONODNGNOD BANONOANONONON BOO NONONONOANO 
—_ ew” wet ue —— — weet eat —_ — —— — — — — — —_ aw wat 

wed met ret et ed ed et LI LINO IN IN ION SS SOMO LM LO LH MOO 
PJ ees os co co J JO am ae Oe) Om OO OL 
bem ee <p eS SE a ee rp 


34 


cS 
LW Fae 
Y <{ 
-— 
= Ls 
eB O 
oe Zz 
a) LW 
ate oc 
OO LW 
ti 
ot LL 
(es) ~~ 
— 6 {LJ 
» © SE 
Cz — 
ae 
— eo YW 
=x< 2a) 
CJ a 
IeV — 
NaS = 
Y= & 
OC ve x< 
we i OC 
ium, 
ZANAIAS 
Orme @¢D 
ZN 
mm Uj LL Lis 
D==z «Ao. 
wines © 
Oaiw tt 
AWZOWMN 
WW LL ee Lo 
OnOdtaze 


DNV ea LOVE 
—a & & ee & eC) 
—eeaeese <j 
MH iloor-e 
—w~ee We ow b= 
bm pe be Ze 
ef id eit 
ee ee) 
(a Mo Mam leis ae 
Soo jG ec@ 
STU Og We (pS Ty ee 
— — 
OarN OT 
SOQ, Oo 


=| ae =! —_— = 





Feed 


",F5.2,2X, REFERENCE ANGLE O 


a 
ry 
on 
amy poet MN 
Tile @ 
ee eM) —— 
Om WI WI ~P 
mt NI bp a a et 
LL wc <f ow ow 
eee ae 
— DSS) rat rent 
_ Nae OOoOw~a 
ace -~er CYrOO 
BT = == WOOrFF 
= 24 MNAOCUWO 
Oo - meet COO OO LI LU 
i Liu cS — >> 
> WZ Z2zIMNww 
< > OW v»srama _ 
SS) — ree LL) peed eg - 
<< iL Pm LS > ee Pr 
<f = m= LEI LU —Yew 
— em) ome fo Sl A a ale Sia ao A a WwY) 
LL - ee te WwW fm LL pa a=; fo @ [eae [pT 
= x «KeZLO WU Le WU OL agitkes 
— mOOOZ tL et rZzw 
Ls CGC Ce CH ned ets LL LL — << 
ed —_ RKPOrFOd =sGO Y 
O <tme <E LS et LL —OQWe2Z2e6 ey BOS OS 
= ZemS>oo>iww Le Ze tOO AWN 
<< ie a SO Wo eee OO << <{ 
OF eee Dev rem i S=szIVNECS ee 
O © ss CD aa) 2a 
Zz Lat LE LU LI LoL | em SSS 
a Eire eEe Er Olt Ze ee SOS 
EN fe eb Ot te Lt Ca a Ye 
- °® 7 Seeeeweeweitiliee were a ie - w= 
om C10) mm OO] Oe em ew ell B&B Bm em Fe eH & Al « & 
NESS LL NICO OK we KOK OK KOK MOK OK OOK OX OK LUI 
el & € BLALALIAIN IAIN KOILALAIANGIOIIAYLY LAAIAIA 
exate [letanie @ ww wewewewee~wMe weewewewe Dw ewe ws @ 
Le LULOOO9O9000 OSB QO00000 AOO00 
Ne Siw Ow ee w= w= @© @ & eereweweeee Dee © 
see ee et et Se ee et ee rt Se eee re et ee Net eet et ee Set Set Set Se 


od mete feed el mee eal eel aed el ced el deed peed ned eed eed eect od md eel eed ed de 
Iaqeadadadtdadddtqadtdae acide 


ZSzrPTPlAATTATASTLTATALTATA TT SzAT TT TA 

ZeMaevweeterrrrertrrrrtertereernrrv~ervoo 
2) (Giese es) ee | eo! Go oe) as | qm | eo te fap) ee | es an | 6 | ow |S [es (a) eb Les | an) 
PPT UL LL tL EL UE LLL LL LLU th LEE LIL LLL iL ud 


= 


NOD ORORAIMNOK DROSUMFOHOHORA 
ES (2 cleo) ololalelelelelelelele(e(slele le lm le, 
eet NINO MAM M QAM MO SEPP STON OOR EE 


55) 


PLANFORM FUNCTIONS 


ETA)) 


= I 
Nie 
|— 
aa 
os! 


rmWM 


on <t | 
aI~O 
—_— e 
LL CO 
T+ 
— 24 
=~ <I oO 
wet LL! ON 
>< OL 


O 
ze 2Q 
Comm 


re{) OZ 
—t OOK 
) aA 2] 
Fl | el 
Se aus 
LL mt ox OK LL 








a ee 


732395*DS INC THE TA) 


36 


SUBROUTINE RMAT(SPL,RP1L»eMESH,AR,RR) 


—m 
sj] 
{WwW 
Or 
eC) pte, 
2a cen — 
fi—— a a) 
Ion Own 
we SII == 
° OF al Bi) — = 
% AN <{<{ 
micah eee )C) 
<{ eCh hy) oad om 
Use OO ALU LU aL LL 
CQO fF NS ZAaANC™ 
COUN mt th HE te 
—- Oe Qe tod 
m=O PR NOON I ot 
OWN! LNA 
mt C9 Ze et NI Et] it 


WWM. OMA Ou 
QE ZN LNNHNIL 
SZ OY) HY} Wk 
mea DQOMNSSFOD 


cS 
he @ 
e #C) 
om vant {{ 
Wot an 
pana oom 


tote 


OOg 
oor 


mit SO. IE 
i oo Tee 
2=7A0 e 

n<ito 
meted OL Ii 

m1) OL 
Oot ioe 
Qoarx 





ae a ee! ia =~ 


HETAP) *XIP 


eZ 
nee <he 
om Slav 
I+ 
—— 
Wa, o 
ran 
| 
~~ IA 
ae) 
Or 
LO 
OWN 
Ix 
ON ee 
ac. il 
pe 
wr = 


PID2+(L+l)/2 


pd 
12 
=o 


L,R 
lyS 
1D2 


lou 


UPA Kk 
Cem 
L=K*S 


oa 


Ta 

_ e am iL) 

i> NI e 

° | re 

oOo = HH + 

=) <a ele 

— <« < os = 

re Oo pS —k FT 

ae eel pee LL , <2 c 

@  t. <r st ~ 

a te Fi 2 

qe cle mMO F#N © 

te) rm <4 =m 

CNS ee OT 

Sorat oF GT wet 
mm OmMEA+— KO Cure 
~WY~eN 1 Owe {ime wwe 3 f= T+ 
ONTtowulea¢eHe 2e wr 
SEO eh RW RW OKm 
ew eet Ne |W RH, % IY) 


set lh we ees KI Te TOOT aSI~O 
hue QOmmR—m 1 3 CL ke Ze 
me ee UO 1] QU er Ly et AI SLU 

Q eG Ot UK TET ow NONH 
MANO-K 1 ICON New| O- 

me <I WH Neen oi 
OC Tim Drm hs me CD) PIA IN 
QOODO8W KE WOTTOTOOVNWNM 


D7 


(DSIN( (K-1)*PHI)*DSIN( PHI )) 


= 
i 
Ce 
VY) 
© 
OO 
OO 
~+ 
INO ~ 
~ @ /) 
UN aad % 
— = 


~— Ile We O 

— mt il 
Lew toOs 
mm TOT Od 


Wm Os 
= 


DTHETA)/(PI*AR) *(RiItR2) 


SZ 
|x 


(PHI-(DSIN(2 .0*PHI}/2.0)) 
*((DSIN( (K-2)*PHI)/(FLOAT(K)I-2.0))-(DSIN(K*PHI )/FLOAT(K))) 


71,8,9 
(PHI+OSINCPHI)) 


x 
PH 


% NS 
— SAIN MUN Ww 
AUD met Qe NI TEL 
es +. sero 
~OoOdDC00O=—V~ 
~NEW ROn 
N NN NZil— 
ieee) ee EIN 
tec) 2 COI oe 


= DO om 
— 


F2 
+C (THETA)-HIP) 


0.5 


10,11,12 


eG 2) 


mt OL ete PE ee ON DZ LUIDD 
OWN) iles SATAN YIN ase NZS H#YAZTZ 
OT ZHaA e ~— <I 3 x LK OO HN wt em eH 
SWUWOO ae Qe STS NANNY YEH RONZERED 
Le SEN I met 1 ott TU UES SRO 
2a——HtOom WODCLZ YOY Ll CY TOnAV~OOowz 
=emOQOWWA GOAret AOY HMTOTONNOOCVOXWw 


east tk 2 ib 

NN N2Zwe2Zz2Zreo 
CSOxvrCtorovodw=a 
COYKOCTOCOOVAU 


QO ms NIN OAS Ay so N= 
— at — =m 


a 
<f 
= 
Lu ae 
i ad pas 
= = 
a (2 
= — 
ae 28 
a's OL 
of =_ Pah 
~~ ~ O — 
—_J ie, — Y) 
Ne <I = a2 OQ 
~ ~~ Nie ff _ N ~ 
= ale OC « LL — * am — — 
< ~ x eeN LT NO onal -~ = 
<< -_ ad -~ TO - el + = ~_= = 
= Lu ce aoe) | = & Oz oe | — (au sae 
te de OQ fone l~z = os 2 
oe N YI wYINA® =%* + + — Cw 
me OO a ~ ol OY Ree) ~w ot J “At gay (2 ok 
CO * (2 di —- On ~~ #H#aAQ Ce Oa +* Iw 
+f AJ + qt ~Aama CS at) = = — mm te if 
NW = = UWA (aa we QO +O — ww # 
YI + (ae So tweet — AHSAN BOW -D | la 
% uy ce Wes IO NN &aOong ! Oo | MO YY ret 
um AN 1 Ww avo NO LNs LwNa “OQ ~~ *##HLI 
‘ e _ a C22 =. eye -A.—-0O XY eQ MNP we UO Il 
GSIinow tutu — te OO As aes ce See QO ZS2iline emit 
= 
= 
© 
fea 
fom) 
aD 
Y) 


58 





SUBROUTINE DGMPRD(A,ByRyNy MeL) 


elololeilojlaleleolelelololelololea 
DOAANMPUOH DAHOSAIG 
COP SP SE NSPS SA PSP LALA LS 


(adja ladle Saas ie Me Me die die Ae Sie ae Ge ae 
eevee oc: a 6. oO. 6.0. 46.0:0- 6.0: 0 4. 
ee ee Pe eo ee ee 
SOOOCOOOL HS OOOOOOOO® 


(am) 
mae 
oO 
% 
oe) 
a) 
—— 
<I 
+ 
42 z= —_ 
» ~ ~ ce 
— od — 


it it Ii 
Ysa = Oma 
+ +2 I +e NZ 
TSFOM~Oe!) x~—mOo~ Oa 
Qf mat tet Ot CHO OD 
wou tf Wouwoum Howe CO 
CYOWw~Ce Owe Oe Owntiz 
ves ee et OD ee OH OOO SO OL 


© 
= 


De) 





dN‘T=I € OO 

YOLISA @ 3HL YOS 3A10S 

(DNS T=IS(CdNS TEP SCF ST IOI (TOESI)DSLIOM 
(OO0C*9)3SLIUM 

(Olav dN) LVWO 11V9 


XIYLVW O SHL JLNdWOD 


WO O 


WOOO 


(dNS T=14(1)dG) 
UV 
HSAWSTdufTds 


we, a 2 oe 
came? “eee Soe SO te 
Ly Uo LL Li LL 
Z RREEF 
de aa de aad 


Vivd 
(dN*T=I5(1)d0) 


ja 
ONO 


eh 


ee ee et ONS et ee 
nl exe | amo U [aaa = By a 


<- 
OO] NNS Fe ANMS 


SO t OOmTo NO D DBOOCSOS 


mew OO LOVIN | He UY 


il 

d 
UV 
dS 


HS AWS Tdd fT 


aM -@°@ <0 <@0pl le a ee ee 


<I 

eS 

~ 
OOUN 


7120050" ¢ *O SOMO. 
(STIIXSISTISONIMS (STIYONIMS (GTIEONIM S(SGT)IZONIMS (G 


peal 2 
ewe, CA] pee om LO\ rm 


O <e YO ewe tre Nhe O 


O Ve YO AROwreINn~wo & 
a) te petit be SP 


QO Wrote 


ZO COX >eZ2ze 
OU NW O Whose D OC 


(S°G)LNIWNS £(S 
((S)T 


(a AT) 


OO _a€artr OO aAexXoO eel 


oO Ht) WY YUnt —_— 
eeaodDS Saar |i 


QO ta YO 2ZBZwwieaemwnedtea BF WUWwWwianw zw 


“, 
) wD WS 


WD) OG OF ANN ° 


O wW WO 
Oo & Oreo 


O 


4409093090 99009910 0909999 999 99909399992999 


Hh 
3 a3 
3 i NOTINGIYLSIG JYNSSAY¥d JSIM 
; INIM WSAGI NV wos ANIT VSGWVD NVIAW SHI 
J 2 


b909999999990005 99 7999009093939 9939000099)02 


OM WwW CON etew 070 


LE = us () Ce QL we tp fh 


OO qo OO sewer IWOOX< ean 
WO FRO YO Se CXKWOOOOAOQR —- CKAKUNNeAe A BEVZsT 


GT 
dg 
Jf 
0 
G 
e) 
g 
S 
2, 


OO Wr OF Ako KY @& 
OO FO O A e 
oS => 


OO 
oy ia4 


6, 


60 








Q,NP,1,1-0E-4,ITER2) 


COMPUTE THE CAMBER LINE 


OUO 


Fo 

pen, 

ce 

C4 

=! 
— 1! 
On "sz 
=m ~ 
eo am 
mam nN 
- a 
a 
Pus Ee 
= Y) 
= ~ 
> —_ 
WY if 
- Zz 
ce om 
<{ ies 
- 2 
a os ~ 
Yn = 
LW od 
= — 
- 2 
7 = 
SS 
oa) 
-~ WY 
= a 
CO wo 
Y) a 


ALL CAMBER( 
RITE (6 302) 
RITE (6,301 


Oz= 


RESULTS 


PLOT 


WOO 


rt od CL 


CL He mt OL eee wee 
CC eam wer HE wee 


eo iN O. 


P(RPL,XIWINGL,ITBL;RIB1L) 


P(RPL XI»WING2,ITBL,RTBL) 
P(RPL,XI,WING3,ITBIL,RTBL) 
P(RP1L,XI,WING4, ITBLyRTBL) 


P(RP1L,XI,WINGS,ITB1L,RTB1) 


aN CSN SNSEMSSESOM= 


oes) 
(a 
| | 
ae 
i ee, a a. 
AAO SIN -—~ 
oF em em & Kale) 
Pt 4 as bt Pg «2 
bt be it 
ae eee mm 
py peg ee eed =ai~es §6@ 
(2 Ei ee a Se 
co AOD Hct 089 4 ct 
HNN” NY) fh. + 
oe | | | | 
}] nam om om om th 1] 


enNOTN USN THe ig 


Ot pet Ges peng Ped feet ~) mt peg CH) met ee |] OY me ome CY ome OE ee OY ome (ae 
CWO NONOANONONSO 


we ere if 


© Sai Pl at wed Hom CA ff es a wo — —" —_ we? eae 


MIODOOOR MG MIMO we I tet et edt edt deed 
222222 ew HO IMO ONION YO UOO 
OC tt ee tes OO Te Oa ee ee et EE 


_ 
a) 


61 


NS 
Mm 


OF CHORDWISE 


NG!) 
“NG. 


-_— & 


URE DISTRIBUTION") 


wd MS 
ee Cle ew LL 
~~ oe & & & & 
om OOS ALS OS OK 
OWN FH ALAN 
ae & @& & & G&S 
—ere wee | ww © 
HKa® i OOas 
—wvepe (/) = wo 
ang? Sat fae wager Maa? “oa 
Ree ZEEE 
ae eta 
2S LTS AS 
(a ge dia UU pe aie aes 
2 | a aay = | oo | tf 
0 He Ws Weg) 9 a 
ol 
Oe MOSw 
Oooo 2800 


ad oe = —_ om 





*,14) 


= o 
© Ww cc 
— ui © 
_ —- 
= ~~ WW 
oa) =a W 
= — > 
ro aa 
—_ wc wm 
Fa) Ce 
_ QO wW 
eS So ft 
= 
= iW Ww 
=~ -_- 2n-~w 
WD * me (> 
mY a OO UL. 
Qw) _ © 
FLW ec wel 
Od re WWoaQ 
<Q <I OW 
We Ze >Y 
S ~f <I 
Wu YF eae 
ee © O 
<_I<{ Wiest Llc 
7 Be eles Ce 
MWY bm be LU 


=> @ cE x<e => & 

oe © emem @f\] & & & 

<x«K NOK < «KX 

LV UY © eNWUn yy 

=> Ox—e es @ee & 
COxXLUOOQOOC ae 
cece =—wOew wf ewe 

—w wey “eae ‘eae “ane! “ang? Seow? See? weet te 

a a a ede ee 
aie daddadadad 

> > i> a i oe 
Seer ceecerwaaOo 
OoCOCoqcoocora 
UL WE LLL UL th ek ULL UL oy 


OM DOH OAINOoS 
COQOoceoeo0osdo 
Se Qi NUM MMOMN Oo 


ON XLCTHETA) 


62 


ION C(THETA) 


UNCT 
=1.0 
ETURN 
ND 


= 
C 
R 
= 


FUNCTION DXL( THETA) 


DXL=0.0 
RETURN 
END 


<I 
j- 
UJ 
ae 
—~ 
<< ee 
- WY 
Li 
oe) 
—_ + 
ww §=CdLSY, 
WwW oo 
a Mm 
GQ 
a ™ 
oO 
—ONZ 
mH © eC 
OOD 
ZulunheG 
wsBOUOWS 
Wien yea e ¢ bo 





SUBROUTINE QMATO(NP, AR, Q) 


_ ° 
i aad ® 
= ur’ NSN 

- e ~ 
tr Oo a ee] 
ra | =O. = 
ONO wD OO + 


mefinwl awe iW OD 
Neate eee Om FOZ 
ae Qo MN oxime 
Woes OHS er HD 
= i — et PTS SRKO 
CACAO TNCOOCCUW 


Nee 


SUBROUTINE CAMBER(SP1,RP1L,MESH,AR,SUMINT »ByNP) 


C\] om ~_ a 
Ow on ee) 
=~ o rw 
Qu Own 
im == 

One ape Te? 
AF = pe be 
OL ee LN << 
NS x eG el | | 
Ser GN 
]NHNMO QULIUILL LL 


GC: ans 
OLN eat Hh OE a 

OOO xtH00 
Cee Pnad il is 
UN Nest il Nerf 
OVZae eNi i De 
WWM) CO = OU 
eet a ae | en © 4 9 ee Be 
Ze OL OL 
SsOOANZzZAO 


63 


) 
(THETAP)*XIP 


xqH% oO Oud 
AON el 
CLNAN je LW 
OL ete we JT 
es~Ot NE 
Bo oe 1 LN 
=ZaAQa eu 

y1<IOn x< 
aoe iO fi 

m™ UO <o 
Oortia}- es 
OO m MLL 


i 
-_ r 
NY & 
e } 
a} a? 
! be 
i ae? 
~~ ee 
—w §«<t 
— OO 
—< Jj _ 
Sue > 4 
J+ -_ + 
LL me 4% _ 
ee a <I 
sl oP) = 
—~mW  w 
moO W”) pS We 
OL me S po 
OW ww — 
=Om—— ©) 
rm ws ( O<T TO + 
rune | Oram 
ONTO SWIe 


ZZO el CWe 
e ele | WW 


Oe unt ane-~_ 


(Oe Met 11 OU ~w 


© il Q ety) ij 
BAMMO t Ore ti ax 
z=. ome HW LL << — I 
DDO ODT es fe We 
NNQOADA*«-FWOTL 





*( SINC CK -LD*EPHI D/C FLOAT (K)-1.0)-SIN((K4+1)*PHI)/(CFLOAT(K) +1 


am LO) aA. 

—™ Ff ll 

l= — iS, = 

=I+ oS t oc 

+ + + ® at <t 
—- dad=— FT N <{ x 
oa~ FN ™ © ~ ~ — 
<I Ut — * ae uu ee 
- Ot = N — eye ~ 
wwi<t FT ~ a OR ~ 
ee st = % GOe < _ 
w TO 2 eS =] ero ke Bate 
Oo A> mm a. 8 mt <u aoa 
Co wN WN -_ MN OLs*e eT =e) 23 > 
et +#+H HX z= ~ Dp LU WI) > 
— amtt ae: ——, Pd _ WAG — Ot Y) 
<x am «J VI x r= ~~ ~ Tt OWY %+ <I + 
— Es ewe YY) we NS mt fe FE oe Ike = 
+ wr w —~C=zO !} o= N wwegqgOunMn KUO 
-~ Tao t-— mitt oe =) = OeskeTleE win 
<q Fest PI LHnwn Lo NY &ND SS <Xws—D Teh 
—- ww whee CAAaAWN+ «& N+ OOND ~OnSN FO 
Wow oll Zt Ht KIiNw He HHEW oH + Twse OXKDW 
Leen ¢€ TT eats + NN HNtee aN 1 aot +N e 
melMOLT< Ree OO—egwew Tea N EO WESwsSTFZzoNudt OF DPWSTVNZTW 
~~ 1% eT eww eet SIM aeNaeN DIDDIDVI=—= ent# eODDD eowD 
ASIHONNWZMY OOYKHW I «© © ZHtnZnnDdD ONnaseonnZ7»0fr2ZzZ 
XT owe MNOeO ti Nts CooOo steht knee ele il tw Hoe 
OnN RT tno ANA ete It Rm Hm NAR AIAN I mee J INR TH LONE 2 
Hew it Way ou == STE ON Nw Sete ee SS | in | Sa i me a 


Mei ettUeste OL DDuU LOT eC DD CD e222 moma OD DD Ow]! 
TOLTTOOWNH VYNMWNODNM He LTOT eOnmMonnnnw LTLIOnNNnOnmnnnocw 


= 
Mm Oo PF N aa) om 
— 


64 





: 


=. 
_ 


REFERENCES 


Glauert, H., THE ELEMENTS OF AEROFOIL AND AIRSCREW 
THEORY, 2nd, ed., Cambridge University Press, 1959. 


Prandtl, L., "Applications of Modern Hydrodynamics 
to Aeronautics", NACA Report 116, 1921. 


Multhopp, H., "Methods for Calculating the Lift 
Distribution of Wings (Subsonic Lifting Surface Theory)", 
bereish Anes, RGMeZooaemel o> O:- 


Perkner, V.%o, “Ihemsoluebon wom ebrt cing bdan Problems. by 
Worcex Lattice [neon erm icmeumc., ROM 7591, 1947. 


Weissinger, J., "The Lift Distribution of Swept Back 
Wings", N&CA Tecnnical Report 1120, 1947. 


merkins,€.E.,Runvyan,H.b., and Woolston,D.S., “On the 
Kernel Function of the Integral Equation Relating in 
Lift and Downwash Distribution of Oscillating Finite 
Wings in Subsonic Flow", NACA Report 1234, 1955. 


Lopez, M.L., and Shen, C.C., "Recent Developments in 
Jet Flap Theory and Its Application to STOL Aerodynamic 
Analysis", AIAA paper number 71-578, revised 19/72. 


Giesing, J.P., Kalman we Go eomemnodden, N.F., “Subsonic 
Unsteady Aerodynamics for General Configurations", 
AIAA paper number 72-26, 1972. 


Kuethe, A.M., and Schetzer, J.D., FOUNDATIONS OF 
AERODYNAMICS, 2nd ed., Wiley, 1959. 


Abbott, I.H., and Von Doenhoff, 4.E., THEORY OF WING 
SeeCUlOns Diver, Fos 7. 


65 





Pe teorSs tT RIBULION ELST 


Defense Documentation Center 
Cameron Station 
Alexandria, Virginia 22314 


i brary, Code 02] 2 
Naval Postgraduate School 
Monterey, California 93940 


Professor T.H. Gawain, Code 57 Gn 
Department of Aeronautics 

Naval Postgraduate School 
Monterey, California 93940 


iedweamo@nm GC. Parks 
COMRESTACSUPWING 

NAS New Orleans 

Belle Chasse, Louisiana 7/0146 


Chairman, Department of Aeronautics 


Naval Postgraduate School 
Monterey, California 93940 


66 


No. copies 
2 











